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Abstract: We present the results of our QCD analysis for polarized quark distribution 
and structure function xgi{x,Q'^). We use very recently experimental data to parameter- 
ize our model. New parameterizations are derived for the quark and gluon distributions 
for the kinematic range x e [10~®,1], e [1,10^] GeV^. The analysis is based on the 
Jacobi polynomials expansion of the polarized structure functions. Our calculations for 
polarized parton distribution functions based on the Jacobi polynomials method are in 
good agreement with the other theoretical models. The values of Aqcd and as{Mz) are 
determined. 
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1. Introduction 

The nature of the short-distance structure of polarized nucleons is one of the central ques- 
tions of present day hadron physics. For more than sixteen years, polarized inclusive deep 
inelastic scattering has been the main source of information on how the individual partons 
in the nucleon are polarized at very short distances. 

The theoretical and experimental status on the spin structure of the nucleon has been 
discussed in great detail in several recent reviews (see, e.g., Refs. ||l|, ^, |^, ^. During the 
recent years several comprehensive analysis of the polarized deep inelastic scattering (DIS) 
data, based on next-to-leading-order quantum chromodynamics, have appeared in Refs. [5- 
29]. In these analysis the polarized parton density functions (PPDFs) are either written 
in terms of the well-known parameterizations of the unpolarized PDFs or parameterized 
independently, and the unknown parameters are determined by fitting the polarized DIS 
data. 

Determination of parton distributions in a nucleon in the framework of quantum chro- 
modynamics (QCD) always involves some model-dependent procedure. Instead of relying 
on mathematical simplicity as a guide, we take a viewpoint in which the physical picture 
of the nucleon structure is emphasized. That is, we consider the model for the nucleon 
which is compatible with the description of the bound state problem in terms of three 
constituent quarks. We adopt the view that these constituent quarks in the scattering 
problems should be regarded as the valence quark clusters rather than point-like objects. 
They have been referred to as valons. In the valon model, the proton consists of two "up" 
and one "down" valons. These valons thus, carry the quantum numbers of the respective 
valence quarks. Hwa and et al. [30-38] found evidence for the valons in the deep inelastic 
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neutrino scattering data, suggested their existence and applied it to a variety of phenom- 
ena. Hwa [^] has also successfully formulated a treatment of the low-px reactions based 
on a structural analysis of the valons. In |4C] Hwa and Yang refined the idea of the valon 
model and extracted new results for the valon distributions. In |41, 43| unpolarized 
PDFs and hadronic structure functions in the NLO approximation were extracted. 

In Ref. the polarized valon model is applied to determine the quark helicity distri- 
butions and polarized proton structure functions in the NLO approximation by using the 
Bernstein polynomial approach. 

The extraction of the quark helicity distributions is one of the main tasks of the semi- 
inclusive deep inelastic scattering (SIDIS) experiments (HERMES El, COMPASS p[. 



SMCllg]) with the polarized beam and target. Recently in Ref. [28| the polarized valon 
model was applied and analyzed the flavor-broken light sea quark helicity distributions with 
the help of a Pauli-blo eking ansatz. The results reported of this paper is good agreement 
with the HERMES experimental data for the quark helicity distributions in the nucleon 
for up, down, and strange quarks from semi-inclusive deep-inelastic scattering |44]. The 
reported results of Ref. |44| is based on the Bernstein polynomial expansion method in the 
polarized valon model framework. 

Since very recently experimental data are available from the HERMES collaboration 
1 47] for the spin structure function gi, therefore there is enough motivation to study and 
utilize the spin structure and quark helicity distributions extracted via the phenomenolog- 
ical model. Since these recently experimental data are just for different value of and 
not fixed Q^, one can not use the the Bernstein polynomial expansion method. Because in 
this method we need to use experimental data for each bin of separately. So it seems 
suitable, to use the Jacobi polynomial expansion method. 



In this paper we use the idea of the polarized valon model to obtain the PPDFs in the 
LO and NLO approximations. The results of the present analysis is based on the Jacobi 
polynomials expansion of the polarized structure function. 

The plan of the paper is to give an introduction to the polarized valon distributions 
in Section 2. Parametrization of parton densities are written down in this section. In Sec- 
tion 3 we present a brief review of the theoretical background of the QCD analysis in two 
loops. The method of the QCD analysis of polarized structure function, based on Jacobi 
polynomials are written down in Section 4. A description of the procedure of the QCD fit 
of gi data are illustrated in Section 5. Section 6 contains a results and conclusions of the 
QCD analysis. 



2. Polarized valon distributions 



The idea of nucleon as a bound state of three quarks was presented for the first time in 



Ref. 1 45]. On the other hand the similar idea, which is called valon model was presented 



in Ref. ]37]. According to the valon model framework, one can assume a simple form for 
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the exclusive valon distribution which facilitate the phenomenological analysis as follows 

GuuD/p{yi,y2, ya) = giymTyl^iyi + y2 + ya - i) , (2.1) 

where yi is the momentum fraction of the i'th valon. The U and D type inclusive valon 
distributions can be obtained by double integration over the unspecified variables 

Gu/p{y) = j dy2 j dy3GuuD/piy,y2,y3) (2.2) 
= gB{a + l,b+l)y%l-yr+'+\ 

Gn/piy) ^ J J dy2GuuD/p{yi,y2,y) (2.3) 

= gS(a+l,a + V(l-y)2°+i. 
The normalization parameter g has been fixed by 

/ Gu/p{y)dy = [ GD/p{y)dy = 1 , (2.4) 
Jo Jo 

and is equal to g = [B{a + l,b + l)B{a + l,a + b + 2)]~^, where i?(m, n) is the Euler-Beta 
function. R.C. Hwa and C. B. Yang |4^] have recalculated the unpolarized valon distribu- 
tion in the proton with a new set of parameters. The new values of a , b are found to be 
a = 1.76 and b = 1.05. 

To describe the quark distribution q{x) in the valon model, one can try to relate the 
polarized quark distribution functions or to the corresponding valon distributions 
and G^. The polarized valon can still have the valence and sea quarks that are polarized 
in various directions, so long as the net polarization is that of the valon. When we have 
only one distribution q{x, Q^) to analyze, it is sensible to use the convolution in the valon 
model to describe the proton structure in terms of the valons. In the case that we have 
two quantities, unpolarized and polarized distributions, there is a choice of which linear 
combination exhibits more physical contents. Therefore, in our calculations we assume 
a linear combination of G'^ and G^ to determine respectively the unpolarized (G) and 
polarized (AG) valon distributions . 

According to unpolarized and polarized valon model framework we have |25] 

qi/p{x,Q') = Y. jyU3il^Q'')Gj/p{y)'^ , (2.5) 

and 

Ag./p(^,Q') =E [' ^Q^/li-,Q')^GJ/piy)^ . (2.6) 

■ J X y y 
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where qi = qj + and Aqi = qj - q\ and also Gj = G] + G] and AGj = G] - G]. 
As we can similarly see for the unpolarized case, the polarized quark distribution can be 
related to a polarized valon distribution in Eq. ( |2.6| ). 

The polarized parton distribution of the polarized constituent quarks in the Eq. 
is dependent on a probe |^]. It is assumed in more detail that the polarized constituent 
quarks distribution, that is, Agj/j(z = in the scale of is equal to A(jfj/j(z, Qg) = 



5{1 — z). From Eq. ( |2.6| ) now is obvious that the polarized quark distribution in the scale 
of Qq is equal to polarized valon distribution. 

Here we want to use the polarized valon distributions which introduced in |25]. Accord- 
ing to improve polarized valon picture the nonsinglet polarized valon distribution functions 
are as following 



AG,.(y) = AW,.(y)xG,(y), (2.7) 

where 

AW,.(y) = i,A,y"^{l - yf^{l + 7,y + , (2.8) 

the subscript j refers to U and D-valons. The motivation for choosing this functional form 
is that the y"^ term controls the low-y behavior valon densities, and (1 — y)^^ terms the 
large-y values. The remaining polynomial factor accounts for the additional medium-y 
values. The normalization constants Aj for U and D-valons are as following 

A^j-^ = [B{a + a^^,,2 + a + b + (3^)+r]^B{Q.b + a + a^,2 + a + b + I3^j) 

+-1^B{1 + a + a^,2 + a + b + p^j)] / B{1 + a,2 + a + h) , (2.9) 

and 

^d"^ = [5(fo + aD,2 + 2a + /3^) +r?^S(0.5 + 6 + a^,2 + 2a + /3^) 

+-f^B{l + b + a^,2 + 2a + l3^)]/B{h+l,2a + 2) . (2.10) 

These quantities are chosen such that the are the first moments of AGj{y), = 
dyAGj{y). Here B{a, h) is the Euler Beta-function being related to the F-function. 
In the present approach the QCD~evolution equations are solved in Mellin-A^ space 
as will describe in section 3. The MELLiN-transform of the valon densities is performed 
and Mellin-A moments are calculated for complex arguments A by y^~^ AGj/p{y)dy. 

As seen from Eq. there are five parameters for each valon distribution. To meet 
both the quality of the present data and the reliability of the fitting program, the number 
of parameters has to be reduced. Assuming SU{3) flavor symmetry and a flavor symmetric 
sea one only has to derive one general polarized sea-quark distribution. The first moments 
of the polarized valon distributions AG^^ and AG^ can be fixed by the SU{3) parameters 
F and D as measured in neutron and hyperon /3-decays according to the relations : 

2C^-C^ =F + D, 

2i^+i^ = ?,F-D . (2.11) 
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The factor 2 in Eq. ( ^.11 ) is due to the existence of two-J7 type valons. A re-evaluation of 



F and D was performed in Ref. [24| on the basis of updated /3-decay constants [49| leading 
to 

= 0.926 ± 0.014 , 

= -0.341 ±0.018 . (2.12) 

This choice reduces the number of parameters to be fitted for each nonsinglet valon density 
to four. 



On the other hand the singlet polarized valon distribution which is defined in |25] is 
as following 

AG'^iy) = AW-iy) X G,{y) (2.13) 
where AWj{y) in Eq. ( |2.13|) has the following form 

5 

AWjiy) = AWjiy) x ^ Amv"^ • (2.14) 

The additional term in the above equation, (^ term), serves to control the behavior of the 
singlet sector at very low-y values and the AVVj(y) defined in Eq. ( |2.8D . 



In the next sections we need to use the Mellin moments of Eqs. ( p.7| , 2.13 ) to deter- 
mine Mellin moments of PDFs in our analysis. 

3. The theoretical background of the QCD analysis 

Let us define the Mellin moments for the polarized structure function g\{x^Q'^): 

g{{N,Q^)= [' x^-'gP,{x,Q^)dx (3.1) 
Jo 

In the QCD-improved quark parton model (QPM), i.e., at leading twist, and to 
leading logarithmic order in the running strong coupling constant as{Q'^) of Quantum- 
Chromodynamics (LO QCD), the deep-inelastic scattering off the nucleon can be inter- 
preted as the incoherent superposition of virtual-photon interactions with quarks of any 
flavor q. By angular momentum conservation, a spin-^ parton can absorb a hard photon 
only when their spin orientations are opposite. The spin structure function gi has then a 



probabilistic interpretation, which for the proton reads [50| 



= i(e2) [AqsiN, Q2) + AqNsiN, Q^)] . (3.2) 

Here, the quantity —Q^ is the squared four-momentum transferred by the virtual photon, 
N is the order of moments, Cq is the charge, in units of the elementary charge |e|, of quarks 
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of flavor q, (e^) = Yliq^'q/^q the average squared charge of the Nq active quark flavors, 
and Aq{N,Q'^) is the quark helicity distribution for quarks of flavor q. Correspondingly, 
Aq{N,Q'^) is anti-quark helicity distributions. Moreover the flavor singlet and flavor non- 
singlet quark helicity distributions are defined as 

AqsiN, Q2) = [A(z(iV, Q') + AqiN, Q^)] 
q 

= AJ:{N,Q^), (3.3) 

and 

Aqr,s{N,Q^)=-\-Y,el [Aq{N,Q^) + AqiN,Q^)] - Aqs{N,Q^). (3.4) 

For the analysis presented in this paper, only the three lightest quark fiavors, q = u,d, s, 
are taken into account and the number of active quark flavors Nq is equal to three. 

The twist~2 contributions to the structure function gi(N,Q'^) can be represented in 
terms of the polarized parton densities and the coefficient functions ACf in the Mellin 
-N space by |] 



9l{N,Q'') = + ^AC^f )[A<z(iV,Q') + Ag-(iV,Q2)] + ^2AC^ Ag{N ,Q'')} , 

(3.5) 

in this equation the NLO running coupling constant is given by 

Inlln^ 



^^°^4= ^ fin. 



MS \ \ 2 

* MS 



where b = ^^^^ and b' = ^^§^^t^- In the above equations, we choose Qq = 1 GeV^ as 
a fixed parameter and A is an unknown parameter which can be obtained by fitting to 
experimental data. 



Aq^{N,Q^) + Aq{N,Q^), Aq{N,Q^) and Ag{N,Q^) are 

7f , AC^ 



In Eq. (U), Aq{N,Q^ 
moments of the polarized parton distributions in a proton. AC^ , AC^ are also the iV-th 



moments of spin-dependent Wilson coefficients given by 



AC. 



and 



N 



-S2iN) + {SiiN)f + ( - 



1 



AC. 



N 



1, 



^ 2' N{N + 1 
with Sk{N) defined as in Ref. 



2 N{N + 1) 
N -1 

■{Si{N) + l) 



,,11 1 
5i(iV) + _ + _ + 



+ 



N'^ 2N N + 1 2 
2 



iV2 N{N + 1)^' 



(3.7) 
(3.8) 
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According to improved polarized valon model framework, determination of the mo- 
ments of parton distributions in a proton can be done strictly through the moments of the 
polarized valon distributions. 

The moments of PPDFs are denoted respectively by: Auy{N, Q^), Ady{N, Q^), AS(iV, Q^) 
and Ag{N, Q^). Therefore, the moments of the polarized u and (i- valence quark in a proton 
are convolutions of two moments: 

Au^iN, Q2) = 2AM^^(iV, Q2) x AM^/p(iV) , (3.9) 

Ad,{N, Q2) = AAf^^(iV, Q2) X AA4/p(iV) , (3.10) 

In the above equation Mj^^{N) is the moment of AGj/p{y) distribution, i.e. AMj^^{N) = 

/Jy^-i AG^/p{y)dy. 

The moment of the polarized singlet distribution (S) is as follows: 

AS(iV,Q2) = AM^(iV,Q2)(2AM;;/p(iV) + AM;/p(iV)) . (3.11) 

here M'.'^^iN) is the moment of AG^yp(y) distribution, i.e. AM'.'^^{N) = y^-^ AG'.i^{y)dy. 
For the gluon distribution we have 

A<7(iV, Q') = AM3^{N, Q^){2AM'(N) + AM'j,,(N)) , (3.12) 



where AM^^{N,Q'^) in Eqs. (|3:9|,|3l0|) , AM^{N,Q'^) in Eq. (|3Tll ) and also AM9i{N,Q 



which is the quark-to- gluon evolution function, are given in Ref. p5| 



It is obvious that the final form for gi{N,Q'^) involves the total of 15 unknown pa- 
rameters. If the parameters can be obtained then the computation of all moments of the 
PPDFs and the polarized structure function (PSF), gi{N,Q'^), are possible. 



4. The method of the QCD analysis of PSF 

The evolution equations allow one to calculate the Q^-dependence of the PPD's provided 
at a certain reference point Qq. These distributions are usually parameterized on the basis 
of plausible theoretical assumptions concerning their behavior near the end points x = 0, 1. 

In the phenomenological investigations of the polarized and unpolarized structure func- 
tions, for example xg^ or xF^ for a given value of Q^, only a limited number of experimental 
points, covering a partial range of values of x, are available. Therefore, one cannot directly 
determine the moments. A method devised to deal with this situation is to take averages 
of the structure function weighted by suitable polynomials. 

The evolution equation can be solved and QCD predictions for PSF obtained with the 
help of various methods. For example we can use the Bernstein polynomial to determine 
PPD's in the NLO approximation to obtain some unknown parameters to parameterize 
PPD's at Qq. In this way, we can compare theoretical predictions with the experimental 
results for the Bernstein averages just in Mellin-A^ space. To obtain these experimental 
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averages from the E143 and SMC data we need to fit xgi{x,Q'^) for each bin in 

separately |25]. 

If we want to take into account very recent experimental data for all range value of 
to determine PPDFs, it is convenience to apply Jacobi polynomials expansion and not 
the Bernstein one. The advantage of application for all data points, especially very recently 
HERMES experimental data |47|, and not just a series of data points, is our motivation to 
use Jacobi polynomials and not Bernstein polynomial to study spin dependent of parton 
distribution function. 

One of the simplest and fastest possibilities in the PSF reconstruction from the QCD 
predictions for its Mellin moments is Jacobi polynomials expansion. The Jacobi polyno- 
mials are especially suited for this purpose since they allow one to factor out an essential 
part of the x-dependence of the SF into the weight function [^] . Thus, given the Jacobi 
moments a„,((5^), a structure function f{x, Q"^) may be reconstructed in a form of the series 
11-11 

N max 

xf{x, Q2) = ^/^(i _ x)" ^ a„(g2)0-,/3(^)^ (4.1) 

n=0 

where N^ax is the number of polynomials and Qn'^{x) are the Jacobi polynomials of order 
n, 

n 

Q':,^P{x) = Y,cf\a,(3)x\ (4.2) 
i=o 

where c^"'' (a, (3) are the coefhcients that expressed through F-functions and satisfy the 
orthogonality relation with the weig ht xP{l — x)'^ as following 

dx x^{l - x)-G"'^(x)e°'^(x) = 5k,i , (4.3) 





For the moments, we note that the dependence is entirely contained in the Jacobi 
moments 



an(Q') = C dxxf{x,Q'')Qt^{x) 
Jo 



n 



^cf (a,/?)/(j + 2,Q2). (4.4) 

j=0 



obtained by inverting Eq. 4T, using Eqs. (4^ , O) and also definition of moments, /(j, Q^) = 
dx x^~'^xf{x, Q^). 

Using Eqs. ( [4.1| - [4.4D now, one can relate the PSF with its Mellin moments 

^m ax ^ 

xfff— (x,g2) = x^(i_x)" ^ G°'^(x)^cf)(a,/3)5i(i + 2,Q2), (4.5) 

n=0 j=0 

where gi{j + 2, Q"^) are the moments determined by Eqs.( ^^j3.5| ). N^ax, o. and (3 have to be 
chosen so as to achieve the fastest convergence of the series on the R.H.S. of Eq. ( [4.5| ) and 
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to reconstruct xgi with the required accuracy. In our analysis we use Nmax = 9, a = 3.0 
and f3 = 0.5. The same method has been apphed to calculate the nonsinglet structure 
function xF^ from their moments [59-63]. 

Obviously the Q^-dependence of the polarized structure function is defined by the 
Q^-dependence of the moments. 

5. The procedure of the QCD fit of data 

The remarkable growth of experimental data on inclusive polarized deep inelastic scatter- 
ing of leptons off nucleons over the last years allows to perform refined QCD analyzes of 
polarized structure functions in order to reveal the spin-dependent partonic structure of 
the nucleon. For the QCD analysis presented in the present paper the following data sets 



are used: the HERMES proton data ||6J, the SMC proton data |£5|], the E143 proton 



data |51], the EMC proton data |6£, 67]. The number of the published data points above 
= 1.0 GeV^ for the different data sets are summarized in Table 1 for data on gi together 
with the X— and Q^-ranges for different experiments. 



Experiment 


x-range 


Q^-range 
[GeV2] 


number of data points 


Ref. 


E143(p) 


0.031 - 0.749 


1.27 - 9.52 


28 


a 


E143(p) 


0.031 - 0.749 


2, 3, 5 (Fixed) 


84 


m 


HERMES (p) 


0.028 - 0.660 


1.01 - 7.36 


19 


a 


HERMES (p) 


0.023 - 0.660 


2.5 (Fixed) 


20 


0] 


HERMES (p) 


0.026 - 0.731 


1.12-14.29 


62 




SMC(p) 


0.005 - 0.480 


1.30 - 58.0 


12 


m 


SMC(p) 


0.005 - 0.480 


10 (Fixed) 


12 


m 


EMC(p) 


0.015 - 0.466 


3.50 - 29.5 


10 


m 


EMC(p) 


0.015 - 0.466 


3.50 - 29.5 


10 


m 


proton 






257 





Table 1: Published data points above Q'^ = 1.0 GeV^. 

In the fitting procedure we started with the 15 parameters selected, i.e. 4 parameters for 
each non-singlet polarized valon distribution, 6 parameters for singlet polarized valon dis- 
tribution and Aqcd to be determined. For this set of parameters the sea-quark distribution 
was assumed to be described according to SU{3) flavor symmetry. 

In the further procedure we fixed 2 parameters at their values obtained in the first 
minimization and chose the first moment of polarized Non-singlet valon distributions in 
the polarized case [^] . The lack of constraining power of the present data on the polarized 
parton densities has to be stressed, however. Since only more precise data can improve the 



situation, so we add the recent experimental data achieved from HERMES group [ 47 1. 

The final minimization was carried out under the above conditions and determined 
the remaining 15 parameters. The values and errors of these parameters along with those 
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parameters fixed in the parametrization, Sec. 2, are summarized in LO and NLO in table 2. 
The results on Aqcd are discussed separately in section 6. The starting scale of the 
evolution was chosen as = 1 GeV^. 

Using the CERN subroutine MINUIT ||6^ , we defined a global for all the experimen- 
tal data points and found an acceptable fit with minimum x^/d.o.f. = 236.577/242 = 0.978 
in the LO case and xV^-O-f. = 225.920/242 = 0.933 in the NLO case. In this table we 
compare the results reported in |2^ which are based on Bernstein approach and the results 
of the present analysis. Also we obtain the uncertainties of the parameters in Jacobi Ap- 
proach which are not calculated in Bernstein approach. We should notice that the Jacobi 
and Bernstein polynomials are merely used as a tool in the fitting procedure and the results 
are independent of it. 





Bernstein Approach [^5| 


Jacobi Approach 


LO 


NLO 


LO 


NLO 




VclIug 




VclIug 


A V ) .l\/r^A/" 


203 


235 


201±110 


245±58 


o^u 
Pv 
7y 
Vu 


0.0020 

-2.3789 

-L7518 

11.0804 

-1.4629 


0.0038 

-2.1501 

-0.8859 

10.6537 

-0.1548 


0.0015(fixed) 
-2.4489±0.024 
-1.9050±0.196 
13.5420±0.378 
-1.8339±0.120 


0.0018(fixed) 
-2.3590±0.028 
-0.9683± 0.206 
22.6095± 0.418 
-1.9778± 0.150 


D \ D ' 

Id 


-0.005 
-1.5465 
-1.8776 
8.5042 
-0.8608 


-0.0046 
-1.5859 
-1.5835 
9.6205 
-0.8410 


-0.0029(fixed) 
-1.6059±0.029 
-2.2009±0.353 
10.2371±0.829 
0.8751±0.091 


-0.0026 (fixed) 
-1.6638±0.034 
-1.6214± 0.453 
13.1966± 0.942 
0.8483±0.142 


Ai 

A2 
A3 

A4 

A5 


0.0004 
0.2954 
-6.9134 
30.9851 
-39.7383 
16.4605 


-0.0025 
-3.1148 
15.8114 
-21.1500 
10.5025 
-0.9162 


0.0003± 10~^ 
0.2611±0.191 
-6.5429±0.174 
29.6999±0.475 
-37.9379±0.202 
15.5879±0.041 


0.0052±10-^ 
-3.0034±0.241 
14.8696± 0.214 
-18.8057± 0.051 
8.4598± 0.298 
-0.3410±0.050 


X^/ndf 


154.98/123=1.260 


115.62/123=0.940 


236.577/242=0.978 


225.920/242=0.933 



Table 2: Parameter values in LO and NLO of the parameter fit for Bernstein and Jacobi 
approaches. 



6. Results and conclusions of the QCD analysis 

We have performed a QCD analysis of the inclusive polarized deep-inelastic charged lepton- 
nucleon scattering data to next-to-leading order and derived parameterizations of polar- 
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ized valon distributions at a starting scale Qq together with the QCD-scale Aqcd in the 
polarized valon model framework. 

The analysis was performed using the Jacobi polynomials-method to determine the 
parameters of the problem in a fit to the data. 

A new aspect in comparison with previous analyzes is that we determine the parton 
densities and the QCD scale in leading and next-to-leading order by using Jacobi polyno- 
mial expansion method. 

Detailed comparisons were performed to the results obtained in other recent parame- 
terizations |2|, |6|, |7|, The previous results are widely compatible with the present 



parameterizations. These distributions can be used in the numerical calculations for po- 
larized high-energy scattering processes at hadron- and ep-colliders. 

Looking at the dependence of the structure function gi{x,Q'^) in intervals of x 
gives insight to the scaling violations in the spin sector. As in the unpolarized case the 
presence of scaling violations are expected to manifest in a slope changing with x. The 
proton data on gi{x, Q^) have been plotted in such a way in Figure 1 and confronted with 
the QCD NLO curves of the present analysis. Corresponding curves of the parameteriza- 
tions [ p4| , |69| , [70| , fn\ , 72 1 are also shown. Slight but non-significant differences between the 



different analyzes are observed in the intervals at low values of x. However, the data are 
well covered within the errors by all analyzes. 

In Figures 2-3 the fitted parton distribution functions, in leading and next-to-leading 
order for all sets of parameterizations 59, 7C, ^ and their errors are presented at the 



starting scale Qq. 

The polarized structure function x(^{x.Q'^) measured in the interval 3.0 GeV^ < < 
5.0 GeV^, Figure 4, using the world asymmetry data is well described by our QCD NLO 
curve. We also compare to corresponding representations of the parameterizations 70 
69, 24 1, which are compatible within the present results. 



In Figures 5-8 the scaling violations of the individual polarized momentum densities 
are depicted in the range x e [10~^, 1], e [1, 10^] GeV choosing the NLO distributions. 
The up-valence distribution xAuy, Figure 5, evolves towards smaller values of x and the 
peak around x ^ 0.25 becomes more flat in the evolution from 

Q2 = 1 GeV^ to Q2 = 

lO^GeV^. Statistically this distribution is constrained best among all others. The down- 
valence distribution xAd^, Figure 6, remains negative in the same range, although it is 
less constraint by the present data than the up-valence density. Also here the evolution is 
towards smaller values of x and structures at larger x flatten out. 

The momentum density of the polarized singlet quark xAS, Figure 7, is positive in 
the kinematic range shown for all and for a; > 0.1 , but changes sign for lower values of 
X. The maximum of the distribution at = 1 GeV^ around x ~ 0.3 moves to x ~ 0.2 at 
= 10^ GeV. At the same time a minimum around x ~ 0.03 moves to x ~ 0.002. 

The momentum density of the polarized gluon xAG, Figure 8, is positive in the kine- 
matic range shown for all Q^. Also in this case the evolution moves the shape towards 
lower values of x and flattens the distribution. By having polarized parton distributions, 
the first moments of the polarized parton distributions can be obtain. The first moments 
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Figure 1: The polarized structure function as function of in intervals of x. The error 
bars shown are the statistical and systematic uncertainties added in quadrature. The data are 
well described by our QCD NLO curves (solid lines) which based on the valon model and Jacobi 
polynomials expansion method. The values of C{x) are given in parentheses. Also shown are 
the QCD NLO curves obtained by AAC (dashed hues) [||, BB (dashed-dotted hues) !§, GRSV 
(long-dashed lines) (t^], LSS (dashed-dashed-dotted lines) @ and DNS (dotted lines) f| for 
comparison. 

of the polarized parton densities in NLO in the MS scheme at = 4GeV^ for different 
sets of recent parton parameterizations are presented in table 3. 



Distribution 



Model 



Ref. 



BB 



GRSV |70| 



AAC m] 



Adv 
Au 
Ad 
Aq 
AG 



0.9207 
-0.3391 
0.8390 
-0.4207 
-0.0817 
0.9325 



0.8769 
-0.3313 
0.8145 
-0.3937 
-0.0624 
0.8322 



0.926 
-0.341 
0.851 
-0.415 
-0.074 
1.026 



0.9206 
-0.3409 
0.8593 
-0.4043 
-0.0625 
0.6828 



0.9278 
-0.3416 
0.8399 
-0.4295 
-0.0879 
0.8076 



Table 3: Comparison of the first moments of the polarized parton densities in NLO in 
the MS scheme at = 4GeV^ for different sets of recent parton parameterizations. The 
second column (Model) contains the first moments which is obtained from the valon model 
and the Jacobi polynomials expansion method. 
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Figure 2: LO polarized parton distributions at the input scale Qq=1.0 GeV^ compared to results 
obtained by BB model (dashed hne) (ISET=1)[|§, A AC (dashed-dotted hne) (ISET=1)[||, GRSV 
(dashed-dotted dotted line) (ISET=3)[|o| and LSS (dashed-dashed dotted line) (ISET=3)[|ll 



We can also obtain the first moment of (the Ellis-Jaffe sum rule) by 

r?(Q2)^ f\xgP{x,Q^). (6.1) 
Jo 

The results have also been given in table 4. The first moments of polarized parton distri- 
butions also shown for some value of Q^. 









AS 


Aq 


A5 


r? 


1 


0.9260 


-0.3410 


0.0965 


-0.0814 


0.5850 


0.1161 


3 


0.9215 


-0.3393 


0.0923 


-0.0816 


0.8651 


0.1195 


5 


0.9202 


-0.3389 


0.0911 


-0.0817 


0.9837 


0.1205 


10 


0.9189 


-0.3384 


0.0898 


-0.0818 


1.1383 


0.1215 



Table 4: The first moments of polarized parton distributions, Au^, Ad^, AS, Aq, Ag and 
in the NLO approximation for some value of Q^. 

In the framework of QCD the spin of the proton can be expressed in terms of the first 
moment of the total quark and gluon helicity distributions and their orbital angular mo- 
mentum, i.e. 

i = IaS^ + A/ + L?, (6.2) 

where is the total orbital angular momentum of all the quarks and gluons. 

The contribution of addition of ^AS and A^ for typical value of = 4 GeV^ is around 

0.978 in our analysis. We can also compare this value in NLO with other recent analysis. 
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Figure 3: NLO polarized parton distributions at the input scale Qq=1.0 GeV^ compared to results 
obtained by BB model (dashed hne) (ISET=3)[||, AAC (dashed-dotted hne) (ISET=3)[||, GRSV 
(dashed-dotted dotted line) (ISET=l)[|o) and LSS (dashed-dashed dotted line) (ISET=l)[|l| 




Figure 4: Polarized proton structure function xgf measured in the interval 3.0 GeV'^ < < 5.0 
GeV^ as a function of x. Also shown are the QCD NLO curves at the same value of obtained by 
LSS (dashed-dashed dotted hne) (ISET=1)@, GRSV (dashed-dotted dotted line) (ISET=1)[0, 
BB model (dashed-dotted line) (ISET=3)p| and AAC (dashed line) (ISET=3)[pl for comparison. 



The reported value from BB model H is 1.096, AAC model IJ is 0.837 and also GRSV 
model @ is 0.785. 
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Figure 5: The Polarized parton distribution xAu^,, evolved up to values of (5^=10,000 GeV^ as 
a function of x in the NLO approximation. The solid line is our model, dashed line is the AAC 
model (ISET=3)(24|, dashed-dotted hne is the BB model (ISET=3)||6| and dashed-dotted-dotted 
line is the GRSV model {ISET=1)||. 

In the QCD analysis we parameterized the strong coupling constant in terms of 
four massless flavors determining Aqcd- The LO and NLO results fitting the data, are 



aJJcd^ = 201 ± llO(stat) MeV, LO, 



a: 



'■QCD 

(4)MS 
QCD 



245 ± 58{stat) MeV, NLO, 



(6.3) 



These results can be expressed in terms of as{M'^): 

as{Ml) = 0.1281 ± 0.0094(siat) MeV, LO, 
a^(M|) = 0.1141 ± 0.0036(stat) MeV, NLO. 



(6.4) 



These values can be compared with results from other QCD analyzes of polarized inclusive 
deep-inelastic scattering data 



E154 |] : 


a,(M|) 


= 0.108 - 0.116 , 


SMC |7|] : 


a,(M|) 


= 0.121 ±0.002(stat) , 


ABFR |1C] : 


a.(M|) 


^-,^^+0.004 , ,+0.009 ,, 
= 0.120 exp) (theor 
-0.005 ^ ' -0.006 ^ 


BB |6S] : 


a,(M|) 


= 0.113 ± 0.004(stat), (ISET=3) 



and with the value of the current world average 



a,(M|) = 0.118 ±0.002 [f4| . 



(6.5) 
(6.6) 



- 15 - 




Figure 6: The Polarized parton distribution xAd^, evolved up to values of Q^ = 10,000 GeV^ as a 
function of x in the NLO approximation. The solid line is our model, dashed line is the AAC model 
(ISET=3)@, dashed-dotted hue is the BB model (ISET=3)|||| and dashed-dotted-dotted line is 
the GRSV model riSKT^n^l 
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Figure 7: The Polarized parton distribution xAE, evolved up to values of (5'^=10,000 GeV^ as a 
function of x in the NLO approximation. The solid line is our model, dashed line is the AAC model 
(ISET=3)[|||, dashed-dotted hue is the BB model (ISET=3)|||| and dashed-dotted-dotted line is 
the GRSV model (ISET=1)0. 



We hope our results of QCD analysis of structure functions in terms of Jacobi polynomials 
could be able to describe more complicated hadron structure functions. We also hope to 
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Figure 8: The Polarized parton distribution xAG, evolved up to values of Q^=10,000 GeV^ as a 
function of x in the NLO approximation. The solid line is our model, dashed line is the AAC model 
(ISET=3)@, dashed-dotted hne is the BB model (ISET=3)g and dashed-dotted-dotted line is 
the GRSV model (ISET=1)01. 



be able to consider the symmetry breaking of polarized sea quarks by using the polarized 
structure function expansion in the Jacobi polynomials. 
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